A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

07 (02.13.4; 08.01.3; 08.16.4; 08.22.3) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



Near-infrared narrow-band photometry of M— giant and 
Mira stars: models meet observations 

R. Alvarez 1 and B. Plez 2 3 

1 GRAAL, Universite Montpellier II, UPRESA 5024/CNRS, F-34095 Montpellier Cedex 05, France 

2 Astronomiska Observatoriet, box 515, S-751 20 Uppsala, Sweden 

3 Atomspektroskopi, Fysiska Institution, box 118, S-221 00 Lund, Sweden, plez@fysik.lu.se 
Received 20 June 1997; Accepted 19 September 1997 



Abstract. From near-infrared, narrow-band photometry 
of 256 oxygen-rich Mira variables we obtain evidence 
about the loops that these stars follow in colour-colour di- 
agrams. We also find a phase lag between indices related to 
molecular band-strength of titanium oxide and vanadium 
oxide. We compute colours for normal M-giants and Miras 
using hydrostatic and hydrodynamic model atmospheres 
and very extensive up-to-date line lists. Normal M-giants 
colours are well reproduced, reaching a high quantitative 
agreement with observations for spectral types earlier than 
M7. The out-of-phase variations of the various spectral 
features of Miras are also acceptably reproduced, despite 
limitations in the modelling. This enables us to confirm 
that the phase lag phenomenon results from the propaga- 
tion of perturbations in the extended atmosphere. It opens 
new perspectives in the spectral modelling of Miras. 

Key words: physical data and processes: molecular data 
- stars: atmospheres - stars: AGB and post-AGB - stars: 
variables: Miras 



1. Introduction 

Red giants and long-period variables provide crucial infor- 
mation on stellar and galactic evolution. It is thus essen- 
tial to be able to interpret their spectra and photometry in 
terms of stellar parameters (effective temperature, chem- 
ical composition, mass, etc) and physical processes (pres- 
ence of shocks, of a wind, etc). In spite of the difficulty to 
model these stars due to the many complications arising 
from low temperatures (molecule and dust formation) and 
the often very extended atmosphere, important progress 
in the modelling of static (Plez et al. 1992; see also the 
review by Plez 1997) and dynamic (Bowen 1988; Fleischer 
et al. 1992; Hofner & Dorfi 1997) cool star atmospheres 



have been achieved in the last years. 
Confrontation with observations is the natural aim of any 
modelling attempt. Checking the validity of the various 
models' assumptions and inputs is obviously essential. 
Some recent investigations have dealt with high resolu- 
tion spectroscopic observations (e.g. Plez et al. 1993) while 
others concentrated on photometry (Bessell et al. 1996). 
They demonstrate that model atmospheres and synthetic 
spectra are now reaching a high quantitative agreement 
with observations, due to the recent great improvement in 
opacity completeness and accuracy. 

The near-infrared narrow-band photometry of Lockwood 
(1972), described in Sect. 2, has various qualities which 
make it very suitable for the study of M-giants and Mira 
variables. It consists of narrow-band filters around 1 /im, a 
spectral region where very cool stars emit most of their en- 
ergy and where, unfortunately, very few observed spectra 
are available. The narrow-band filters used allow a precise 
study of absorption features of molecules dominating the 
optical spectrum of cool stars, titanium oxide (TiO) and 
vanadium oxide (VO), for which new data have been re- 
cently published. Furthermore, Lockwood's data represent 
one of the very few set of regular observations on whole 
cycles of a large number of Miras. The study of the ef- 
fects of pulsation on absorption features is thus possible. 
Observational evidence for a phase lag between indices re- 
lated to TiO and VO absortion features is given (Sect. 3). 
In Sect. 4, synthetic colours from static models are com- 
puted and compared to M-giant colours. Thanks to the 
updated millions of spectral lines taken into account, a 
remarkable agreement is obtained. In Sect. 5 we compute 
LTE radiative transfer in hydrodynamic models that sat- 
isfactorily reproduce the colour variations and the phase 
lag for Mira variables. 



2. The data 



Send offprint requests to: alvarez@graal.univ-montp2.fr 



Alvarez & Plez: NIR photometry of M-giants and Miras 



2 

2.1. Lock-wood's photometric system 

Lockwood (1972; hereafter L72) observed 292 M- and S- 
type Mira variables from 1969 to 1971 at the Kitt Peak Na- 
tional Observatory, with a five-colour narrow-band photo- 
metric system based on Wing's 27 colours system (Wing 
1967). He obtained 1795 individual sets of five-colour mea- 
surements, at different phases, and, for some stars, during 
several cycles. The narrow-bands were chosen to measure 
depths of molecular band heads of TiO and VO, or re- 
gions relatively free of molecular absorption. Table 1 lists 
the properties of the five-colour system. The 78 and 88 fil- 
ters were used to measure several bands of TiO. The VO 
absorption was measured by the 105 filter. The 87 and 
104 filters were expected to match 'continuum' regions. 
However, it appeared to Lockwood that the 87 filter be- 
comes contaminated by TiO bands in the later M stars. 
In L72 are tabulated the measurements of the 104 magni- 
tudes and of the four colours: 78—87, 88—87, 87—104, and 
105—104. The median standard errors were 0.012 mag for 
104 and 0.006 mag for the colours. 



Table 1. Properties of Lockwood's (1972) five-colour system 


Filter 


Peak 


Half-power 


Feature 


designation 


wavelength 


bandwidth 


(Lockwood 1972) 


78 


7818 A 


90 A 


TiO 


87 


8777 A 


82 A 


Continuum (+TiO) 


88 


8884 A 


114 A 


TiO 


104 


10351 A 


125 A 


Continuum 


105 


10506 A 


loo A 


VO 



Due to the regularity of the observations, this set of 
data allows a precise study of the behaviour of TiO and 
VO absorption with phase. For homogeneity reasons, only 
the oxygen-rich Miras have been subsequently considered. 
They represent a sample of 256 stars for which 1501 sets of 
five-colour measurements were obtained at several phases. 
The period distribution of the sample corresponds per- 
fectly to the one obtained for the M-Miras listed in the 
General Catalogue of Variable Stars (GCVS; Kholopov 
1985, 1987). The stars were observed preferably near max- 
imum and minimum of light curves. 

2.2. Molecular band-strength indices 

In order to study the behaviour of TiO and VO absorption 
features with phase, we have considered four indices: 

78—88: this index is proportional to the ratio between two 
different band heads of TiO. It is a good effective tem- 
perature indicator for temperatures above 3000 K and 



thus was prefered to 88—87 originally given in L72. For 
lower temperatures, the TiO bands begin to saturate. 

78—87: the 87 and 88 filters overlap slightly. As a conse- 
quence, the 78—88 and 78—87 indices are very similar. 
The absolute differences are always less than 0.2 mag. 

87—104: Lockwood originally hoped to obtain a colour 
temperature determination from this index as both fil- 
ters were thought to measure regions free of strong 
molecular absorption. But filter 87 proves to be con- 
taminated at least by TiO bands. 

105—104: it is related to VO bands. This index depends 
strongly on effective temperature below 3000 K where 
the VO molecular bands begin to appear. 

The first three indices can be considered to measure 
mainly TiO features (although we will show below that 
VO bands are also present in these passbands) and 
105—104 may be seen as a pure VO index. The indices 
78—88 and 105—104 complement each other very well as 
temperature indicators. They have been used to determine 
effective temperature for a large number of Miras (Alvarez 
& Mennessier 1997). The four colours are each defined so 
as to increase linearly with decreasing effective tempera- 
ture during a cycle. 

3. Variations of TiO and VO absorption features 
during a cycle. Observational evidence for phase 
lag 

3.1. Phase-colour diagrams 

The phases of all sets of measurements are given in L72. 
They were determined from the mean periods given in 
the GCVS and the AAVSO annual predictions. Phase 
zero marks the maximum of the visual light curve. Fig- 
ure 1 shows the four indices calculated from the entire 
data (1501 measurements) plotted against the phase. The 
curves represent the running average computed for 100 
points on a l/16th cycle interval. The data in the phase in- 
terval [0.5 ; 1.0] are reproduced in the interval [-0.5 ; 0.0]. 
Arrows indicate the extremes of each curve. 

Such a plot has of course no real physical significance 
since the data are composed of observations of 256 dif- 
ferent Miras, which exhibit a large spread of periods and 
mean spectral types. Nevertheless, it gives interesting in- 
sights into the global behaviour of the indices with phase. 
As expected, the colours reach a minimum around phase 
zero when the effective temperature is highest and reach 
a maximum around phase 0.5, at minimum light. 
The important feature is the phase lag that exists between 
the 105—104 index and the other indices. The VO index 
has its minimum value at phase <p = 0.1 or so, whereas 
the indices 78-88, 78-87, 87-104 have their minimum 
slightly before phase zero. The same phase shift of about 
0.1 is also found between the maxima: it occurs at tp = 0.6 
for 105—104, and around ip = 0.5 for the other colours. 
Note that the dispersion around the mean value is mini- 
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Fig. 1. The four colour indices as a function of phase (1501 measurements for 256 O-rich Miras). The curves are the running 
average values and the arrows indicate the extreme values 



mum near phase zero for 87—104 and 105—104 and max- 
imum near phase 0.5, whereas the opposite is true of the 
two other indices. This is explained by the sensitivity, or 
the lack thereof, of the various indices in the effective tem- 
perature range above or below 3000 K (see Fig. 7). As a 
consequence, the positions of the maxima of the TiO in- 
dices (78—87, 78—88) which saturate around 1.8, are not 
easy to locate accurately. The VO index (105—104) has 
values clearly greater than zero only near the minimum 
light. Thus, the presence or not of a phase lag is not to- 
tally obvious in Fig. 1. Colour-colour diagrams enable us 
to confirm that feature. 

3.2. Colour-colour diagrams 

Figures 2 and 3 display the average colour-colour varia- 
tion with phase for 78-88/105-104, and 78-87/87-104 
respectively. The average curves of Fig. 1 are simply re- 
ported in those colour-colour diagrams. The filled circles 
mark phase zero and the arrows indicate the direction of 
variation during the cycle. 



In the 78-88/105-104 diagram (Fig. 2), a loop is clearly 
described during the cycle. To a single 78—88 value 
two different 105—104 values corresponds, depending on 
whether the Mira is rising or declining. This loop confirms 
that a phase lag exists between 78—88 and 105—104. On 
the other hand, there is no loop in the 78—87/87—104 di- 
agram (Fig. 3) . It demonstrates the absence of phase shift 
and thus the occurence of extremes at more or less the 
same phase, as can be seen in Fig. 1. Similarly there is 
no phase lag between 78-88 and 78-87. Only 105-104 is 
out-of-phase by an amount of 0.1. 

The presence (absence) of loop in the 78—88/105—104 
(78—87/87—104) diagram is also verified when only a part 
of the sample of Miras is considered. Early and late spec- 
tral type Miras or short and long period sub-samples have 
been tried. The apparent phase shift and the presence of 
a loop in a colour-colour diagram which is its counterpart 
are independent of period and spectral type. The loop for 
shorter periods (higher T c g) is at smaller colour indices 
than that for longer periods (lower T e s). 
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Fig. 2. 78—88/105—104 diagram. The curve represent the run- 
ning average values shown in Fig. 1. The filled circle marks the 
phase zero and the arrow indicates the direction of variation 
during the cycle 
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Fig. 3. Same as Fig. 2 for 78-87 and 87-104. 

3. 3. Individual phase shifts 
3.3.1. Method 

The results described above concern only mean values. 
The trends revealed by this global approach should be ver- 
ified by considering the individual sets of observations. In 
the following, we adopted the complementary individual 
approach to determine the possible phase shifts. 
For each Mira with a number of observations at least equal 
to six (93 objects), we fitted the four indices as a function 
of phase <p by sine curves: 



a = a Q + a x sm[27r(v3 - ip )] 



(1) 



where c, is one of the four colours; ao, a-i and ipo are the 
three free parameters. 

We only used the Miras which have observed colours dis- 



tributed in a phase range larger than 0.3 in order to prop- 
erly fit the observations. We excluded the stars for which: 



max(ci)fi t > max (q) observed + 0.4 



min(c;)fit < min(Q) observed - 0.4 



(2) 



(3) 



For these stars, the extremes given by the fit are too dis- 
tant from the observations and thus are less reliable. Fi- 
nally, we rejected the Miras with: 



max(105 - 104) fit < 0.05 



(4) 



as it means that these stars show no appearance of VO 
absorption bands during the whole cycle (the temperature 
is always too high) and the fits have no meaning in these 
cases. There is no such problem with the other colours. 
Among the 93 Miras with a number of observations greater 
than 5, 76 finally remained. The phase shift between two 
indices is simply: 



Sip = (f (ci) - <A)( c i) 



3.3.2. Results 



(5) 



Figure 4 presents the histogram of phase shifts between 
the (78—88, 105—104) indices. A majority of Miras clearly 
displays a phase lag. The distribution is almost symmet- 
rical, the mean value being 0.08, with a sigma of 0.07. 
The variation of the 105—104 index lags behind that of 
78—88, as it is observed in Fig. 1. The individual phase 
lag values are not correlated with period or spectral type. 
The histogram for the (78—87, 87—104) indices is shown 
in Fig. 5. No global phase lag is apparent as in Fig. 4. The 
calculated mean value is 0.01 with a sigma of 0.07. 
The fit of the observations by sine curves is of course not 
really accurate and the individual phase lags determined 
with this method are not reliable, but both distributions 
provide an excellent statistical confirmation of what was 
concluded from the previous global approach: the 105—104 
index appears out-of-phase with the other indices by an 
amount of almost 0.1, which corresponds to more or less 
30 days for a typical Mira. 

3.4- Comments on the phase lag 

Lockwood (1972) already noted that Miras perform large 
loops in a colour/spectral- type diagram. He pointed out 
that the mean phase of maximum light at 1.04 /im occurs 
at 0.11, whereas at V it occurs at phase zero by definition. 
Mean minimum light at 1.04 /im is reached almost simul- 
taneously with mean minimum visual light, although the 
individual phases of these events often differ by 0.1 cycle 
or more. 

We have shown that the phase lag and the loop phe- 
nomenon are also displayed by the particular near infrared 
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shock fronts in the atmosphere creates a complex tem- 
perature, density, and velocity stratification. Parts of the 
atmosphere are heated up while others cool down. The run 
of temperature is often far from the smooth decrease found 
in static LTE models. As a consequence, during a fraction 
of the cycle, the VO molecular lines can still be building 
up strength, as e.g. VO molecules are forming (the VO in- 
dex increases), while the TiO lines have started to weaken 
(TiO indices decrease): a loop then appears in the corre- 
sponding colour-colour diagram. We will try to investigate 
this hypothesis in Sect. 5 by using a set of dynamical at- 
mosphere models to compute the colours variation during 
a cycle, but this requires that we are first able to well 
reproduce Lockwood's colours for non-variable M-giants. 
This is the aim of the following section. 



Fig. 4. 

93o (78- 
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Fig. 5. Histogram of phase shifts: opo(87 — 104) minus 
<M78-87) 



indices considered here. More precisely, the phase shift is 
between the VO index on one hand and indices mainly 
related to TiO on the other hand. The phase lag between 
changes in the various bands which are dominated by dif- 
ferent molecular absorption (TiO and VO) may be ex- 
plained by the fact that these bands are formed at various 
depths in the atmosphere. What we see is a consequence 
of the shocks (created by the pulsation) running through 
the atmosphere. 

Important progress in nonlinear stellar pulsation models 
has been made in the last years (Feuchtinger & Dorfi 1994, 
1996a, 1996b; Bessell et al. 1996; Ya'ari & Tuchman 1996). 
Dynamical atmosphere models have also been successfully 
developed recently (e.g. Bowen 1988; Fleischer et al. 1992; 
Hofner & Dorfi 1997): the effects of pulsation generated in 
the interiors are simulated by periodic mechanical driving 
at the inner boundary (piston model). All these numer- 
ical simulations show that the propagation of successive 



The first step in any attempt to investigate the origin 
of the phase lag is to correctly reproduce Lockwood's 
colours for M-giants. It is now possible thanks to impor- 
tant progress in the modelling of red giant atmospheres 
achieved in the last years (Plez et al. 1992, hereafter 
PBN92; see also the review by Plez 1997). 

4-1. The synthetic spectrum program 

We computed synthetic spectra with TurboSpectrum, an 
enhanced version of the Spectrum package developed at 
Uppsala observatory and used in previous investigations of 
AGB stars (e.g. Plez et al. 1993). Input model atmospheres 
were calculated with SOSMARCS (PBN92), one of the 
latest developments in the MARCS suite of programs ini- 
tiated by Gustafsson et al. (1975). SOSMARCS generates 
spherically symmetric, hydrostatic, flux constant (convec- 
tive + radiative) models in LTE. It is especially designed 
for cool star atmospheric conditions. It features a full 
Opacity Sampling treatment of all opacities with 10912 
sampling points between 990 and 125000 A. TurboSpec- 
trum is based for a large part on the same ensemble of 
routines (e.g. evaluation of continuous opacities and solu- 
tion of the transfer equations). Any kind of temperature- 
pressure stratification may be used and the execution time 
is relatively short even when millions of lines are taken into 
account. It offers the possibility for instance to test com- 
plex profiles which can mimic pulsating atmospheres, and 
perform extensive calculations with various abundances 
and/or line lists. 

4-2. The atomic and molecular data 

We included for the spectrum computations atomic line 
data from N to Ni (Kurucz 1992). We also included line 
lists for the TiO a, f3, 7, 7', <5, e, <j> and a-f electronic 
transition systems, the VO A-X and B-X electronic tran- 
sition systems, and a newly calculated line list for H 2 
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(J0rgcnsen, private communication). The absorption sys- 
tems of 12 C 14 N and 13 C 14 N (J0rgensen & Larsson 1990) 
are also taken into account, with a 12 C/ 13 C ratio equal 
to 19, which is a typical value after the first dredge- up 
(Smith 1990). 

Most of the TiO and VO band strengths were updated 
since PBN92. In PBN92 the electronic transition momenta 
(i? 2 ) of Davis et al. (1986) were adopted for the TiO a, 
/?, 7, 7', 5 and <fi systems, supplemented by Brett's astro- 
physical calibration for the TiO e (Brett 1990). In subse- 
quent works (e.g. Plez et al. 1993) the lifetime measure- 
ments of Doverstal & Weijnitz (1992) replaced some of the 
earlier values. Recently Hedgecock and coworkers (1995) 
performed lifetime measurements with a better accuracy 
and for a larger number of electronic states than previ- 
ously for the a, (3, 7 and 7' systems. These new lifetime 
values are significantly different from the older ones. We 
decided to adopt them in our investigation. We generated 
also a line list for the a-f system using the lifetime given 
by Hedgecock et al. (1995, experimental), and Schamps et 
al. (1992, theoretical) which agree within 5%. This system 
absorbs mostly in the V band region and has no impact 
on the present calculations. For the TiO cf) (Hedgecock et 
al. did not measure the corresponding lifetime) and the 
e (the same authors provide only an upper limit of the 
strength) systems, the results of the recent ab initio cal- 
culations performed by Langhoff (1997) are adopted. His 
calculations for the other bands are in good agreement 
with the Hedgecock et al. values. We included all 5 sta- 
ble isotopes of Ti with terrestrial abundance ratios in the 
computation of the line list. 

The VO A-X and B-X band strengths in PBN92 result 
from the astrophysical calibration by Brett (1990) and are 
not expected to be very reliable especially as they were 
calibrated using outdated mean opacity models. Recently, 
Karlsson et al. (1997) measured the lifetimes of the A, 
B and C electronic states. We included their values in 
our line list which was also revised using Cheung et al. 
(1982a, 1982b, 1994), and Adam et al. (1995). The new 
line positions are shifted by up to 65 cm -1 (—40 A) in the 
B-X system, and 25 cm -1 in the A-X system compared 
to PBN92. The line-strengths are scaled by a factor 0.51, 
0.28, and 0.72 for the A-X, B-X, respective C-X transi- 
tions relative to PBN92's values. 

We summarize the i? 2 and the electronic oscillator 
strength (,f -values) we used in Table 2 and 3 for the 
TiO and VO systems respectively. Like Larsson (1983), 
we prefer to avoid the discussion of / c -values which lack 
a sound and unique definition in molecular spectra. We 
provide them anyway as a guide for comparison with 
other works (e.g. J0rgensen 1994). We define / c through: 
/ e = 3.038 x 10~ 6 x croo x R 2 e x 2 ' whcrc ° m is 

some average wavenumber for the (0-0) band of the tran- 
sition system, / and // denote the upper and lower state 
respectively, and A is the projection of the electron orbital 



angular momentum (0 for £ states, 1 for n states etc); see 
Larsson for more details. 



Table 2. Adopted electronic transition momenta for the TiO 
systems 



System 


p2 

n e, Davis 


R 2 

^c,this work 


resulting 
fe- value 


R 2 

c,this work 

B 2 . 
c, Davis 


a 


3.0 


1.8 a 


0.106 


0.6 


P 


5.14 


2.3" 


0.125 


0.45 


7 


3.6 


1.84" 


0.0786 


0.51 


7' 


2.8 


1.9" 


0.0935 


0.68 


5 


1.4 


1.4 


0.048 


1.0 


t 




0.064 b 


0.0023 




<P 


0.9 


0.32 b 


0.0178 


0.36 


a-f 




1.69" 


0.098 





a based on lifetime measurements (Hedgecock et al. 1995) 
b based on lifetime calculations (Langhoff 1997) 



Table 3. Adopted electronic transition momenta for the VO 
systems 



System 


r>2 a 
n e, Brett 


r>2 b 
c,this work 


A-X 


0.207 


0.106 


B-X 


2.590 


0.735 



a astrophysical calibration by Brett (1990) 

b based on lifetime measurements (Karlsson et al. 1997) 

As shown in Table 2 and 3, almost all TiO and VO 
band strengths have been obtained from recent laboratory 
measurements or calculations. Apart from the 5 system, 
all bands have been updated since PBN92. The choice to 
keep the Davis et al. value for the S system comes from 
comparisons between observed and calculated spectra for 
a sample of M giants and dwarfs. Adopting Langhoff's 
lifetime for b 1 !! and branching ratios for the S and <f> tran- 
sitions leads to obviously too strong S bands. We could 
not find an explanation for this mismatch and decided 
therefore to keep the seemingly better value from Davis 
et al. (1986), which also matches Ramsey's (1981) mea- 
surements of the (5(0-0) bandhead at 8160 A rather well 
down to Toff ~ 3400 K. Adopting Langhoff's value, which 
is twice as large, results in a calculated Ramsey's Dsseo 
index systematically too large by 0.1 to 0.14. This is how- 
ever puzzling and will be further investigated by Plez et 
al. (in preparation) in their discussion of the new large 
grid of Uppsala models. 

Davis et al. (1986) used the lifetime of the upper state 
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of a transition in the j3 system measured by Feinberg & 
Davis (1977) to put their measurement of relative transi- 
tion momenta on an absolute scale. Our revision of the (3 
system strength, following Hcdgecock et al. lifetime mea- 
surement, amounts to a downward scaling by a factor 0.45 
of Davis et al. value. For the other bands quoted by Davis 
et al. the factor is between 0.36 and 0.68 (see Table 2). 

4-3. Influence of the various molecular bands on colours 

We first computed individual spectra for each absorption 
band system in order to establish its particular influence 
on Lockwood's passbands. Table 4 presents the systems 
which are really dominant in each band and those for 
which the influence is less pronounced but still measur- 
able. The systems that are not listed have no significant 
influence. 

The most prominent features measured by the filters are 
the TiO 7 and 8 systems, and the VO A-X and B-X sys- 
tems. A noticeable result is the relatively important influ- 
ence of VO absorption on the 87 and 88 filters. Further- 
more, when the T e g reaches values lower than 2800 K, the 
features due to the VO B-X system become even stronger 
than those due to TiO bands in our models. This effect 
was not expected by Lockwood. Concerning the contin- 
uum points, the 104 filter can effectively be considered as 
a good measurement of a region relatively free of strong 
molecular absorption; this is certainly not the case for the 
87 filter. The CN and H2O bands and the atomic lines 
have no significant effect on the emergent fluxes in the fil- 
ter passbands. 

The Lockwood colours, as they measure the influence of 
a limited number of molecular band systems as shown in 
Table 4, constitute an ideal opportunity to check the cal- 
ibration of some of those TiO and VO bands and/or the 
accuracy of the models and the synthetic spectra. 



Table 4. Dominant molecular features in each band of Lock- 
wood's photometric system 



Filter 


Dominant systems 


Less influent systems 


78 


TiO 7 


VO B-X, TiO 7', 8 and e 


87 


TiO 7 and VO B-X 


TiO e 


88 


TiO 8, 7 and VO B-X 


TiO e and VO A-X 


104 




TiO <f>, 8, e and 7 


105 


VO A-X 


TiO e, 7, <f> 



4-4- Comparison to observations 

Lockwood has observed 61 standard red giants with spec- 
tral types from K5 to M8 in the same five-colour system 



as the Mira variables. In order to compare these observa- 
tions with synthetic colours, we used fifteen SOSMARCS 
models representing a sequence of red giants with solar 
composition, a stellar mass of 1.5 Mq and T ff ranging 
from 2500 to 3900 K. The surface-gravities logg vary from 
1.60 at T cff = 3900 K to -0.50 at T cS = 2500 K. These 
models were used by Fluks et al. (1994) and are further 
described there. The emergent fluxes were computed be- 
tween 7500 and 11000 A with a resolution of 0.2 A. Chang- 
ing the sampling to 0.8 A changes the colours by at most a 
few hundredth of a magnitude. Figure 6 shows an example 
of a synthetic spectrum (T e g = 3500 K) . The dotted line 
is part of an optical spectrum of HD 123657 (BY Boo), a 
M5 -giant, obtained by Serote Roos et al. (1996) at a reso- 
lution of 1.25 A using the Aurelie spectrograph, attached 
to the OHP 1.52m telescope. The region covered ends at 
8920 A. The agreement is very good. 
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Fig. 6. Synthetic spectrum with T cff = 3500 K, log #=0.9, 
M=1.5 Mq. The dotted line is part of an optical spectrum of 
HD 123657, a M5-giant, from Serote Roos et al. (1996). The 
Call triplet lines are not included in the synthetic spectrum. 
The filter transmission curves used in this work are superim- 
posed (full line) 

Lockwood's filters of Table 1 were applied to the fifteen 
spectra. The zero point calibration is based on a spectrum 
of Vega (Dreiling & Bell 1980) as used by Bessell & Brett 
(1988). An effective temperature has been assigned to each 
standard giant using the spectral types given in L72, and 
the spectral type-T e ff relation determined by Fluks et al. 
(1994). The overall good agreement of the synthetic and 
observed colours is demonstrated in Fig. 7. 

A small but systematic discrepancy exists for the 
78—88 and 78—87 indices when they are greater than 0.8 
(Fig 7a andb). Furthermore, both colours are not ever- 
increasing along the sequence of T e g. The observed vari- 
ation with phase of these indices in the Miras sample 
brings the conclusion that this behaviour is apparently 
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Fig. 7. Comparison of synthetic and observed indices in temperature-colour diagrams. Crosses are the red giants observed by 
Lockwood. Diamonds represent the sequence of hydrostatic models with T e g ranging from 2500 to 3900 K 



not real. The synthetic 87—104 colours are maybe also 
slightly above the observed ones for temperatures lower 
than 3300 K (Fig 7c). The 105-104 index appears to be 
the least well- modelled (Fig 7d), but note that its scale of 
variation is a factor 10 smaller than the other colours. 
Discrepancies might be due to uncertainties in the: 

— narrow-band filter shapes (a gaussian shape was as- 
sumed) and in the zero point calibration: this may im- 
ply a systematic error. It seems this is the case at least 
for the 105 filter (see Fig 7d). Standard errors quoted 
by Lockwood for 104 are ±0.012 mag. Errors in our 
zero points are probably of the same order of magni- 
tude. This is problematic only for 105—104, which has 
a very small range of variation. 

— spectral types attributed by Lockwood: some stars 
have the same spectral type whereas they exhibit 
rather distinct colours. Furthermore, the scale of spec- 
tral subtypes (MO to M10) adopted by Fluks et al. for 
their SpT-T e ff relation might differ from the one used 
by Lockwood. 



— laboratory measurements and calculations of lifetime 
as quoted by the different authors: the uncertainties 
are still important despite great progress made re- 
cently in the TiO and VO bands calibration. 

Furthermore, red giants naturally exhibit a spread of 
effective temperature, but also of stellar masses, grav- 
ities and chemical composition. Figure 7 shows only a 
'one-dimensional' fit (in T e g): M. and log<? have a non- 
negligible effect on colours and part of the discrepancy 
might also be explained by this. Test calculations using 
models taken from the grid of Plez (1992) show that two 
models with the same T e g and luminosity but with dis- 
tinct log g and M have differences of tenths of mag for 
some indices (e.g. at 3000 K, the \ogg-— 0.7 and Ai = l 
model differs by -0.06 in 78-87, -0.007 in 78-88, 0.37 
in 87-104, and 0.042 in 105-104 from the logg^O.OO and 
A4=5 model). We recall that the sequence of models was 
computed using the Ridgway et al. (1980) temperature 
calibration, and a theoretical red-giant branch from Lat- 
tanzio (1991) to relate the effective temperature and lu- 
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minosity, assuming a mass of 1.5 Mq. 
The effect of relative abundance changes in C and N 
was also checked. Some spectra were recomputed with a 
change of —0.2 dex in the C abundance and +0.5 dcx in 
the N abundance, which are typical values after dredge-up 
(Smith 1990). The difference is less than a tenth of a mag- 
nitude in all colours at 3100K, and largest for 87—104. 
The fit for the 105—104 index is not excellent below T c g ~ 
3100 K. However, the variation scale is less than 0.2 mag. 
Uncertainties in the observations are of the order of 0.006 
mag for the colours, according to Lockwood. We suspect 
the flux level of our models around 1/im to be slightly 
too high. This appears also when comparing them to ob- 
served spectra of M-giants. There may be a missing opac- 
ity in our computations. This will be further investigated 
by Plez et al. (in preparation) . Note also that the points at 
T cff < 3100 K consist of X Her, EP Aqr, BK Vir, RT Vir, 
SW Vir and RX Boo, which are all variables. In conclu- 
sion, keeping in mind the above discussion, the agreement 
between the synthetic and observed colours is remarkable 
in Fig. 7, and provides a quite safe basis for the study of 
Miras. 

5. Miras and hydrodynamical model colours 

Now that we are confident about the line lists and our 
synthetic colour calculations we may try to compute them 
for Mira stars. An examination of Fig. 1 readily shows 
that the variation of the 87—104 and 105—104 colours for 
Miras greatly exceeds the range of static models, whereas 
this is not the case for 78—87 and 78—88. Together with 
the phase lag between 105—104 and the other colours, this 
shows that static models will not be able to reproduce the 
colour variations along the cycle of Mira stars. 

S. Hofner (1997) kindly provided us with a sequence 
of models (50 in total distributed on 2 successive cycles) 
with the parameters: Al* =1 Mq, L* = 7000 L Q , T* = 
2880 K, P = 390 d, and a piston velocity of 4 km.s -1 . 
This model was already used by Loidl ct al. (1997) to 
synthetize C-rich IR spectra and is similar to the mod- 
els presented by Hofner & Dorfi (1997) except for the gas 
opacity which was adjusted by comparison with standard 
model atmospheres to give more realistic gas densities. 
It was computed with C/0=1.8, which is not adapted 
to our O-rich objects. However, there arc no suitable O- 
rich hydrodynamical models available, and we only want 
to qualitatively understand the phase lag. A likely expla- 
nation is that this phase lag arises from the difference in 
formation depth of the various spectral regions dominated 
by continuum, TiO or VO bands, combined with the run- 
ning temperature T- and density p- perturbations along 
a cycle. For the purpose of this work we may assume that 
the chemistry used in the computation of the hydrody- 
namic model is not a key issue (especially as it is hardly 
reflected in the radiative transfer, treated in the grey case 
with a constant opacity of 10~ 2 cm 2 .g _1 ; the C-rich chem- 



istry has however a strong impact on the dust formation). 
We therefore only used the T and p stratification of the 
models in R, recomputed the chemical equilibrium with 
a solar composition, to derive the electronic pressure P c 
and an optical depth scale r. We did not use all the 500 
points of the original models but selected about 65 depth- 
points with an average separation of 0.2 in A logr, which 
is sufficient for treating the radiative transfer. We slightly 
extrapolated the models inwards to increase the maximum 
t value in order to provide a better inner boundary condi- 
tion (diffusion approximation), but the exact stratification 
deep in the model has no strong impact on the emergent 
spectrum. 

As noted above, we surmise that the phase lag is due to 
a difference in the depth of formation of VO bands (dom- 
inating at 105) and TiO bands. A first and quick quali- 
tative check is offered by the variation in column density 
of TiO and VO along a pulsation cycle of the model. The 
column densities are obtained by calculating the chem- 
ical equilibrium at each depth point of the dynamical 
models and integrating the number densities radially. The 
column densities oscillate up and down along the cycle 
(Fig. 8). The phase lag is 0.065 at minimum column den- 
sity and 0.030 at maximum, the VO variation lagging be- 
hind TiO's, providing a first qualitative confirmation of 
our hypothesis. The ratio of the column densities is about 
50, a factor 5 larger than the abundance ratio of Ti and 
V. 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

phase 



Fig. 8. TiO and VO column densities as a function of phase, 
computed assuming chemical equilibrium in the hydrodynami- 
cal model of Hofner. Both column density variations have been 
scaled to the [ ; 1 ] interval to facilitate comparison. The col- 
umn densities vary between 3.9 10 17 and 1.6 10 18 g.cm for 
TiO, and 5.2 10 15 and 5.0 10 16 g.cm" 2 for VO 

We also computed LTE synthetic spectra for these 
models in the same fashion as for static models. The ve- 
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locity fields were neglected. The resulting colours are dis- 
played in Fig. 9. To facilitate the comparison with observa- 
tions we have also plotted in this figure the mean relation 
from Fig. 1 for each colour. We have allowed for 2 shifts 
of the observed colours, relative to the observations. The 
first shift, identical for all 4 colours, is in phase (the cal- 
culated curves are shifted by 0.08). There are no a priori 
reasons for the zero of the phase to be the same in the 
observations and in the models. For the observations, we 
recall that the zero is defined as the maximum light in V 
band. In the models, the zero is defined as the time when 
the piston moves outwards with its maximum velocity. We 
also shifted the data vertically to allow an easier compar- 
ison of the shape and amplitude of the curves. The scale 
on the right-hand side of the diagrams is for the observed 
data, the one on the left-hand side is for the synthetic 
colours. The agreement is surprisingly good. 

The minima are all reasonably well matched with the 
possible exception of the 87—104 index. From the calcu- 
lated minima we find a phase lag of 0.09 of 105—104 rel- 
ative to 78—87 and 78—88 in very good agreement with 
the value derived from observations (0.08). We also find 
a slight phase lag for 87—104 of 0.04, which is also de- 
tectable in the observations at minimum light, and due to 
the hybrid nature of this index. The maxima are more dis- 
crepant, especially for 78—87 and 78—88, with a tendency 
for the synthetic colours to show a secondary minimum. 
We suspect this to be due to the same cause as the S-shape 
shown in Fig. 7 by the synthetic colours of static models 
around 3000 K; we believe this S-shape to be in error. It 
is possible that additional absorbers missing in our syn- 
thetic spectra would solve the discrepancy. A more serious 
problem is that the 78—87 and 78—88 synthetic indices are 
1.5 to 2 mag too large in absolute value, although their 
amplitude of variation is about right. As 87—104 is of the 
right magnitude, it could be that this problem stems from 
an inadequate value of the 78 flux. A more detailed dis- 
cussion is not useful given the approximations we have 
made. We have used one specific C-rich hydrodynamical 
model stratification for which we have computed an LTE 
radiation transfer with an O-rich equilibrium chemistry. 
We have compared the result of these calculations for a 
unique set of model parameters to mean colour variations 
resulting from the average of 256 different Mira stars. It is 
actually amazing that we get such a good agreement. But 
one should certainly perform the same computations for 
a series of models with various M±, L+, P, M, and check 
for possible variations of the phase lag, of the amplitude 
of variations of the various indices, etc. It would also be 
desirable to use O-rich models. One more limitation of 
our approach is that we have neglected the velocity field 
that shifts spectral lines differentially with depth. In this 
particular model the velocity of the gas is of the order of 
+ 15 to -15 km.s -1 in the line formation zone. This will 
spread out the opacity, block more flux between the lines 



and allow more photons to escape within the lines. This 
may solve some of the discrepancies we discussed above. 

6. Conclusions 

Using Lockwood's (1972) five-colours observations of 256 
O-rich Miras we found evidence for the loops described by 
these stars in colour-colour diagrams, and a phase lag be- 
tween the VO and TiO indices, reflecting the shocks run- 
ning through the extended stellar atmosphere. We used 
the best current hydrostatic and hydrodynamic model at- 
mospheres, and line lists to compute synthetic colours 
for normal M-giants and Miras. We found that the nor- 
mal M-giant colours are rather well reproduced, at least 
for T e ff above 3100 K. Missing opacity sources, non-LTE, 
non-thermal heating, variability may all explain the worse 
match to observations at later spectral-types. This is 
worth further study, but for earlier spectral-types model 
atmospheres and synthetic spectra are reaching a high 
quantitative agreement with observations. This is very en- 
couraging for the next (and last) generation of classical 
models in preparation in Uppsala that will include even 
more up-to-date and more complete opacities (Plez et al. 
in preparation). 

Our experiments in solving LTE radiative transfer in 
one series of hydrodynamical models (Hofner, private com- 
munication) have worked beyond our expectations. We 
found approximately the right variation of the colours and 
the right phase lag between them. The fit is of course not 
perfect, and should not be. We used one specific model 
and the observations are an average of 256 different stars. 
There are many more limitations in the modelling (e.g. 
LTE, no account of velocity field, ...), but the main point is 
that hydrodynamical models, equilibrium chemistry, and 
LTE radiative transfer seem to reproduce at least some of 
the observed data on Miras (cf. also the work of Loidl et al. 
1997 on IR spectra of C-rich Miras). This may be surpris- 
ing at first sight. We cannot just assume that equilibrium 
and LTE prevail all through Miras atmospheres. This has 
to be checked, and that is a heavy task. The advantage of 
the present approach is that it allows us to account for mil- 
lions of spectral lines, carry a large number of tests (vary 
chemical composition, ...) while keeping the computation 
time into reasonable bounds. We could easily compute 
center-to-limb intensity profiles as well. We believe this 
work opens the prospect of calculating accurate colours 
and spectra for cool stars, not least long-period variables 
(LPV's). The out-of-phase variations of the various spec- 
tral features, reflecting the propagation of perturbations 
(shocks) in the atmospheres are a very strong constraint 
to the models. The confrontation to spectrophotometric 
observations, resolving the absorption bands, along the 
whole cycle, will be necessary to test the models and better 
understand LPV's atmospheres, pulsation and mass-loss. 
It will also unveil limitations/inadequacies in the hypothe- 
ses made in the modelling. This would be a very useful 
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Fig. 9. Comparison of synthetic colours of dynamic models (crosses; left-hand side and bottom scales) and observed mean 
relations from Fig. 1 (dotted curves; right-hand side and top scales) 



program for a small telescope equiped with an R « 10 000 
spectrograph, ideally reaching wavelengths up to 1.1 mi- 
cron. 
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